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We investigate the observational constraints on the oscillating scalar field model using data from 
type la supernovae, cosmic microwave background anisotropies, and baryon acoustic oscillations. 
According to a Fourier analysis, the galaxy number count TV from redshift z data indicates that 
galaxies have preferred periodic redshift spacings. We fix the mass of the scalar field as m^> = 
3.2 x 10 _31 /i eV such that the scalar field model can account for the redshift spacings, and we 
constrain the other basic parameters by comparing the model with accurate observational data. We 
obtain the following constraints: ST2 m , = 0.28±0.03 (95% C.L.), Q^o < 0.035 (95% C.L.), £ > -158 
(95% C.L.) (in the range £ < 0). The best fit values of the energy density parameter of the scalar 
field and the coupling constant are fi^o — 0.01 and £ = —25, respectively. The value of Sl^o is 
close to but not equal to 0. Hence, in the scalar field model, the amplitude of the galaxy number 
count cannot be large. However, because the best fit values of fi^o an d £ are not 0, the scalar field 
model has the possibility of accounting for the periodic structure in the N-z relation of galaxies. 
The variation of the effective gravitational constant in the scalar field model is not inconsistent with 
the bound from observation. 



I. INTRODUCTION 

A widespread idea in cosmology is that the universe 
is homogeneous and isotropic above a certain scale. 
This hypothesis, usually called the cosmological princi- 
ple (e.g., [l[), is thought to be a generalization of the 
Copernican principle that "the Earth is not in a central, 
specially favored position" . The assumption is that any 
observer at any place at the same epoch would see essen- 
tially the same picture of the large scale distribution of 
galaxies in the universe. 

However, according to a Fourier analysis by Hartnett 
& Hirano 0, the galaxy number count N from redshift 
z data (N-z relation) indicates that galaxies have pre- 
ferred periodic redshift spacings of Az = 0.0102, 0.0246, 
and 0.0448 in the Sloan Digital Sky Survey (SDSS) (see 
Fig. 2 of 3]), with very similar results from the 2dF 
Galaxy Redshift Survey (2dF GRS) (see Fig. 17 of H) 
Q. These redshift spacings have been confirmed by 
mass density fluctuations, the power spectrum P(z), and 
Npairs calculations ||. The combined results from both 
surveys give characteristic periods of 31. 7± 1.8 h~ x Mpc, 
73.4 ± 5.8 h~ x Mpc, and 127 ± 21 h~ x Mpc @. That 
is, the redshift space for relatively high galaxy number- 
count and other that exhibits comparatively low number 
counts appear alternately. 

127 hr 1 Mpc is the same scale as that found in a pencil- 
beam survey of field galaxies p|. Furthermore, the pe- 
riodicity as a function of z in the distribution of QSO 
spectra has also been reported [1, @] • 

A natural interpretation is that concentric spherical 
shells of higher galaxy number densities surround us, 



with their individual centers situated at our location. 
However, if this interpretation reflected the actual phys- 
ical concentration of galaxies existing at certain dis- 
tances from us, it would definitely be incompatible with 
the cosmological principle that presumes uniformity and 
isotropy of our space-time. In fact, it has been demon- 
strated @, from many numerical simulations using the 
Einstein-de Sitter and ACDM models, that the proba- 
bility of getting such a periodic spatial structure from 
clustering and cosmic web filaments is less than 10 -3 . 

One viable explanation for the periodicity of redshift 
space (the picket-fence structure of the N-z relation) 
that preserves the cosmological principle would be in- 
troducing a scalar field that non-minimally couples with 
the curvature scalar, as originally proposed by Morikawa 
(1990 ,1991) [|, El and extended later by Hirano et al. 
(2008) 11]. Such coupling can produce a periodic struc- 
ture of redshift space in the following way. These scalar 
field models [9i-fll1] assume a potential V(4>) proportional 
to 4> 2 so that the scalar field oscillates around the po- 
tential minimum, inducing an oscillation of the Hubble 
parameter, due to the curvature coupling of <f>. An epoch 
of relatively rapid expansion and an epoch of relatively 
slow expansion appear alternately due to this oscillation. 
During the former epoch, the number density of galaxies 
diminishes noticeably, whereas during the latter epoch, 
the decrease in the number density is less drastic. This 
would then give a periodic structure as an apparent or 
illusionary effect in the N-z relation, because we are ob- 
serving the two types of epoch alternately with increasing 
z. Since the spatial distribution of galaxies remains ho- 
mogeneous for a given instance of time, the cosmological 
principle can still be preserved. 
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In scalar field models |9h11|]. the scale of the periodic 
redshift spacings is determined by the mass parameter 
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m<f, of the scalar field. In this paper, we first fix the mass 
of the scalar field as = 3.2 x 10 _31 /i eV such that the 
scalar field model can account for the redshift spacings of 
the \27h~ 1 Mpc period @j. Next, we constrain the other 
basic parameters of the scalar field model of Hirano et 
al. (2008) II] by comparing the model with accurate 
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observational data from type la supernovae (SN la) 12 1, 
cosmic microwave background (CMB) anisotropies 13], 
and baryon acoustic oscillations (BAO) [lij]. Finally, we 
test the validity of this model on variation of the effective 
gravitational constant. 

This paper is organized as follows. In the next section 
we summarize the evolution equations for the scalar field 
and the scale factor of the oscillating scalar field model 
of Hirano et al. (2008) [ll[ and describe the stationary 
state of the scalar field. In Section ITO! we investigate the 
observational constraints on the scalar field model [ll[ 
using data from SN la [L2\ . CMB anisotropies [HI, and 
BAO [3]. For this model, we investigate variation of 
the effective gravitational constant in order to test the 
validity of the model. Finally, the conclusions are given 
in Section [TV] 



II. OSCILLATING SCALAR FIELD MODEL 

A. Basic Equations 

The action S for our cosmological model is given by 
111: 



where the symbol "□" signifies the d'Alembertian, the 
symbol "(• • -) ;M " is the covariant derivative of the quan- 
tity "• • •" with respect to the spatial coordinate x^, and 
the quantities T^™ r ^ are the (/i, i/)-components of the 
energy-momentum tensor of matter and radiation. 

Note that the potential term (the fourth term on the 
right-hand side of Eq. (TT])) reflects our modification [ll[ 
to the original form employed by Morikawa (1990, 1991) 
[9lHoj. which coincides with the latter if q = 0. This new 
form of the potential enables us to control the epoch when 
the growth of the scalar field ip begins to take place. The 
closer this epoch is to the present era, the less affected 
are the amplitudes of the spatial power spectrum of the 
CMB temperature anisotropy in the large-scale domain 
because of the reduced effect of the late-time integrated 
Sachs- Wolfe effect |ll|. In the observational constraints 
of this paper, because we want to know the value of the 
basic parameters of the scalar field model, we fix the 
constant parameter as q — 0. 

By varying the action (Eq. (TT])) with respect to ■0, we 
obtain the time evolution equation for the scalar field: 

4> = -2^ - 6f ^ - a 2 (^f) 2 (1 - # 2 )0exp (-# 2 ), 
a a \ n J 

where we have employed a dimensionless quantity <p in- 
stead of tp: 
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where ip is the scalar field, R is the scalar curvature, £ is 
the coupling constant, is the mass of the scalar field, 
q is a scalar constant, g^ v are the (p, z/)-components of 
the metric tensor, G is Newton's gravitational constant, 
c is the velocity of light, A is Einstein's cosmological con- 
stant, and L( mr ) is the Lagrangian due to matter and 
radiation. 

Based on the principle of least action, we extremizc 
the action integral by setting 6S — 0, which gives the 
following gravitational field equations: 
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We also employ the short-hand notation A 
where t c is the conformal time defined as 
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dA/dt c , 



f c . , : dA a dA 

t c = -dt, so that A = — = - — , (5) 
J a at c c dt 

where t is the proper time. 

Assuming the Robertson- Walker line element, from 
the (0,0)-components and (l,l)-components of Eq. ((2]), 
and using Eq. ©, we obtain the following set of time 
evolution equations for the cosmic scale factor a: 



66W>+{(6^) 2 + (l-6^ 2 ) 
x ((0) 2 + a 2 (^)% 2 exp(-# 2 ) 
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where i/o is the Hubble constant. The quantities r2 m ,o, 
f2 r o, and £Ia,o are the present mass densities of mat- 
ter, radiation, and the cosmological constant A nor- 
malized, respectively, to the present value of the crit- 
ical density p c ,o(= 3H§/8ttG), viz., O m . = p m ,a/p c ,o, 
r2 r7 o = Prfi/Pcfi and f2yi,o = (c/Ho) 2 A/3. The matter 
included in our model universe is assumed to consist of 
baryonic matter as well as cold dark matter (CDM). 

The corresponding density p^ and the pressure for 
the scalar field <f> are given by diagT^ = (p, p, p, p) with 
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From Eqs. ([8} and ((9|) the scalar field satisfies the 
continuity equation: 
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There is no mutual energy exchange between the mat- 
ter, radiation, and scalar fields. The scalar field interacts 
with gravity merely by coupling with the scalar curva- 
ture. 

The density parameter of the scalar field (f> for the 
present epoch is defined as Q^,o = P<f>,o/Pc,o- Further- 
more, since we are concerned only with a flat geometry, 
we obtain the following constraint from Eq. (J6j) for the 
present era: 



B. Numerical Computation 

Salgado et al. (1996) [H] and Quevedo et al. (1997) 
[lU found that the magnitude of the scalar field must re- 
main extremely small during the epochs around the Big 
Bang nucleosynthesis in order not to effect the success- 
ful prediction made for the cosmic abundances of light 
elements based on the Einstein-de Sitter model universe. 
Further, according to Hirano et al. (2006) (23|, unless the 
value of <j> is kept stationary at almost for a sufficiently 
long time in the early stage of expansion, the WMAP 
observation of the CMB angular spectrum cannot be re- 
produced because of the enhanced effect of the late-time 
integrated Sachs- Wolfe effect. Therefore, in what fol- 
lows, we consider only models having a stationary state 
of the scalar field <j> in the early epoch of the universe. In 
the case of having a stationary state of the scalar field, 
because the value of the scalar field (j> is kept stationary 
at almost in the early epoch of the universe, the ef- 
fective gravitational constant G e ff in the era of the big 
bang nucleosynthesis is almost equal to the value of the 
Newton's gravitational constant G. 

To create such a stationary state of <j> at those periods 
of time without any fine-tuning of the present-day value 
oi_d>o such as those carried out by Salgado et al. (1996) 
[2l| , we integrate the evolution equations in the direction 
from past to present (natural direction). 

A sample behavior of <j> for the oscillating scalar field 
model of Hirano et al. (2008) [TTJ] is shown in Fig. [T] as 
a function of scale factor a. We have used O m ,o = 0.28, 
^.o = 0.01, £ = -25, m = 3.2 x 10~ 31 /i eV, q = 0, and 
H = 72 km s _1 Mpc _1 . 
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Measurements of distant type la supernovae (SN la) 
[l5l IT(| indicate late-time accelerated expansion of the 
universe. For the mechanism for the acceleration, there 
are currently many models including dark energy [13, Gil 
or modified gravity 0, ■ In our scalar field model, it 
is always the cosmological constant which supplies nearly 
all of the driving force to the acceleration. 



FIG. 1. An example of the variation of the scalar field <j> as 
a function of scale factor a, starting from a stationary state. 
The parameters are £l m ,o = 0.28, f2^,o = 0.01, £ = —25, 
m = 3.2 x 10~ 31 /i eV, q = 0, and H = 72 km s" 1 Mpc" 1 . 

Fig. [^indicates the behavior of the Hubble parameter 
normalized to the Hubble constant H/Hq as a function 
of a. The parameters are the same as these of Fig. [TJ 
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FIG. 2. Oscillation of the Hubble parameter normalized to 
the Hubble constant H/Ho. Top right panel is a close-up 
view. The abscissa is the scale factor a. The parameters are 
the same as these of Fig. [T] 



Fig. [3] shows the effective gravitational constant G e g 
given by [| 
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where G is Newton's gravitational constant. The abscissa 
is the scale factor a. 
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FIG. 3. Oscillation of the effective gravitational constant nor- 
malized to Newton's gravitational constant G e g/G. The ab- 
scissa is the scale factor a. The parameters are the same as 
these of Fig. Q] 

In the early epoch, the scalar field starts from a sta- 
tionary state. Then, it begins to oscillate (Fig. Q]). The 
Hubble parameter and the effective gravitational con- 
stant also begin to oscillate almost simultaneously (Fig. 
[H Fig. [3]). The amplitude of the oscillation of the scalar 
field has been damped until the present time (Fig. [T]). 



III. OBSERVATIONAL TESTS 

In this section, we study the cosmological constraints 
on the scalar field model of Hirano et al. (2008) flll j ob- 
tained from SN la [12, CMB anisotropics [11], and BAO 
[l4| observations. Finally, for this model, we investigate 
variation of the effective gravitational constant in order 
to test the validity of the model. 

In scalar field models ll], the scale of the periodic 
redshift spacings is determined by the mass parameter 
to^ of the scalar field. We fix the mass of the scalar field 
as m<f, = 3.2 x 10~ 31 h eV, such that the scalar field model 
can account for the redshift spacings of the 127 7i _1 Mpc 
period [2; . Because we want to know the value of the 
basic parameters of the scalar field model, we fix the 
constant parameter q introduced by Hirano ct al (2008) 
[HI as q = 0. The case of q = coincides with the model 
by Morikawa (1990, 1991) [1,113, Kashino & Kawabata 
(1994) [Hi, Fukuyama et al. (1997) [Hj], and Hirano et 
al. (2006) [H. We fix the Hubble constant as H = 72 
km s^Mpc -1 [27j]. 

For the other parameters of the scalar field model [ll| , 
we carry out a detailed investigation of the allowed pa- 
rameter region using the following observational data. 



A. Observational Data 

1. Type la Supernovae 

We first use the SN la constitution dataset, which in- 
cludes 397 SN la |H]. The 90 SN la from the CfA3 
sample with low redshifts are added to the 307 SN la 
union sample [28j . 

The dataset gives the distance modulus at redshift 

Hobs • 

The theoretical distance modulus (for a given model) 
is defined by 



fi(z) = 51og 10 D L +/i , 



(13) 



where Dl is the Hubble free luminosity distance given 
by 
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and /iq is 
Mo = 5 log 



510 (f^) 



25 = 42.38-51og 10 /i, (15) 



where h is the Hubble constant Ho in units of 
100 km s" 1 Mpc^ 1 . 

We minimize the statistical % 2 function (which deter- 
mines the likelihood function of the parameters) of the 
model parameters. For the SN la data we have 
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where is the total uncertainty of the distance modulus 
0. 

Because the nuisance parameter (Eq. (|15p ) is 
model- independent, we analytically marginalize it as fol- 
lows: 



where 



and 
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We use the XgNj a function in combination with that of 
the CMB and BAO data. 



2. Cosmic Microwave Background 



where the parameter A is calculated by 
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and zbao = 0.35. The % 2 function for BAO is 
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To determine the best value and the allowed region of 
the parameters, we use the maximum likelihood method 
and minimize the following quantity: 



2 _ 2 i 2 i 2 

X — XsNia + Xcmb + Xbao- 



Numerical Results 



(26) 



We now present our main results for the constraints 
from the observational data described in the previous 
subsection. 

In Fig. HI we plot the probability distribution of the 
energy density parameter of matter f2 TOj o for the oscillat- 
ing scalar field model [HI from the combination of the 
SN la, CMB, and BAO data, where the parameters fi^o 
and £ are marginalized. 



The CMB shift parameter is one of the least model- 
dependent parameters extracted from the CMB data. 
Because this parameter involves a large rcdshift behavior 
(z ~ 1000), it gives a complementary bound to the SN 
la data (z < 2). The shift parameter R is defined as 



H(z) 



dz, 



(21) 



where zqmb is the redshift at recombination. Although 
zcmb depends on the matter density £l m o and on the 
baryon density Sl^o at the ~ 2 percent level, we fix this 
redshift to the 5-year WMAP maximum-likelihood value 
of z CMB = 1090 [3. 

Using the 5-year WMAP data of R obs = 1.710 ± 0.019 
[3, the x 2 function for CMB is 



Xcmb 



R- 1.710 
0.019 



(22) 
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FIG. 4. ID probability distribution of the energy density 
parameter of matter f2 m ,o for the oscillating scalar field model 
from the combination of the SN la, CMB, and BAO data. 



3. Baryon Acoustic Oscillation 

Observations of large-scale galaxy clustering provide 
the signatures of the BAO. We use the measurement of 
the BAO peak in the distribution of luminous red galaxies 
(LRGs) observed in the Sloan Digital Sky Survey (SDSS) 
(l4| . which gives 



A obs = 0.469 x 



0.96 
098 



-0.35 



±0.017, 



(23) 



The best fit value is il m .o = 0.28. This is similar to 
the value indicated in the ACDM model. We obtained 
the comparatively stringent constraint as follows. 



n, 



0.28 ±0.03 



(95% C.L.) 



(27) 



Fig. [S] shows the probability distribution of the en- 
ergy density parameter of the scalar field fi,* o f° r the 
oscillating scalar field model [ll| from the combination 
of SN la, CMB, and BAO data, where fl m ,o and £ are 
marginalized. 
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FIG. 5. ID probability distribution of the energy density FIG. 6. ID probability distribution of the coupling constant 

parameter of the scalar field O^o for the oscillating scalar £ for the oscillating scalar field model from the combination 

field model from the combination of SN la, CMB, and BAO of the SN la, CMB, and BAO data, 
data. 



The best fit value of the energy density parameter of 
the scalar field fl^o = 0.01, but not 0. 

We obtained the very stringent constraint as follows. 



ft^o < 0.035 



(95% C.L. 



(28) 



We find that the energy density parameter of the scalar 
field fi^.o must be small in the universe at present. 
Fig. [B] shows the probability distribution of the cou- 



pling constant £ for the oscillating scalar field model [11 1 



from the combination of the SN la, CMB, and BAO 
data, where Q m ,o and are marginalized. Note that 
Fukuyama et al. (1997) [26| pointed out that the use 
of negative values for the coupling constant £ yields a 
succession of mini-inflationary (concave shaped) cosmic 
expansions with time, thereby prolonging the cosmolog- 
ical age. We therefore employ negative values for £. We 
search by step of logarithm for parameter £. 

The best fit value of the coupling constant is £ = —25. 
The constraint is 

£ > -158 (95% C.L.) (in the range £ < 0) (29) 

The absolute value of £ cannot be very large. 

Fig. Fig. 13 and Fig. M show the 2D probabil- 
ity contours for the oscillating scalar field model (llj . 
The dotted (green) and solid (red) contours show the la 
(68%) and 2a (95%) confidence limits, respectively, from 
a combined analysis of the SN la, CMB, and BAO data. 
The other parameters are marginalized. 

The best fit value of fl m ,o is 0.28, and the best fit value 
of O^o is close to but not equal to 0. The absolute value 
of £ cannot be very large. 

The amplitude of N from z data (the picket-fence 
structure of the N-z relation) mainly depends on fi^o 
and £ in the oscillating scalar field model. The observa- 
tional amplitude of the galaxy number count has a high 
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FIG. 7. 2D probability contours in the (f2 m ,o, O.^^) -plane for 
the oscillating scalar field model. The dotted (green) and solid 
(red) contours show the la (68%) and 2a (95%) confidence 
limits, respectively, from a combined analysis of the SN la, 
CMB, and BAO data. 



uncertainty, but in the scalar field model, since the most 
likely value of £l<p.o is small, the amplitude cannot be 
large. 

However, the scale of the periodic redshift spacings 
is determined by the mass parameter of the scalar 
field. The scalar field model can account for the redshift 
spacings for the 127/i _1 Mpc period @, assuming a mass 
of the scalar field of = 3.2 x 10~ 31 h eV. Because the 
most likely values of Q^o and £ are not 0, the scalar field 
model has the possibility of accounting for the periodic 
structure in the N-z relation of galaxies. 

Fig. [TU] shows the normalized energy density of radia- 
tion tt r , matter f2 m , cosmological constant J7a, and scalar 
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FIG. 8. 2D probability contours in the (Q m ,o, £)-plane for the 
oscillating scalar field model. The dotted (green) and solid 
(red) contours show the la (68%) and 2a (95%) confidence 
limits, respectively, from a combined analysis of the SN la, 
CMB, and BAO data. 



FIG. 10. Normalized energy density of radiation Q r , matter 
fl m , cosmological constant Qa, and scalar field Qd>, versus 
the redshift z in the oscillating scalar field model [HI ] with 
the best fit parameters: Q, m ,o ~ 0.28, fl^o = 0.01, £ = —25. 




0.04 0.06 

FIG. 9. 2D probability contours in the (fl^o, £)-plane for the 
oscillating scalar field model. The dotted (green) and solid 
(red) contours show the la (68%) and 2a (95%) confidence 
limits, respectively, from a combined analysis of the SN la, 
CMB, and BAO data. 



field fli/y versus the redshift z in the oscillating scalar field 
model 11] with the best fit parameters: O m ,o = 0.28, 



n 



a 



0.01, £ = -25, where fl 4 



ItiG/ZH 2 



the energy density of scalar field defined by Eq. ([8j . We 
find that the normalized energy density of scalar field f2^ 
is comparatively small in all era. 

In Table U we list the best fit parameters, the x 2 values 
(Eq. (|26))). and the differences of the Akaike information 
criterion (AIC) 29] and the Bayesian information crite- 
rion (BIC) [H, for the scalar field model and the ACDM 
model, from a combined analysis of the SN la, CMB, and 



BAO data. The definitions of AIC and BIC are 

AIC = —2 lni + 2k, (30) 



BIC 



-21nL + /cln N, 



(31) 



where L is the maximum likelihood, k is the number 
of free model parameters, and N is the number of data 
points used in the fit. The % 2 value for the scalar field 
model is smaller than that for the ACDM model. How- 
ever, the values of AIC and BIC for the ACDM model 
are smaller than that for the scalar field model, because 
the number of free parameters of the scalar field model 
is two more than that of the ACDM model. 



C. Variation of the Effective Gravitational 
Constant 

In order to test the validity of the model, we compute 
the variation of the effective gravitational constant. 

From Eq. (| 12[) . the variation of the effective gravita- 
tional constant is 



G 



12£#' 



G cff 1 - 6^ 2 



(32) 



where prime designates the derivative with respect to the 
proper time t. 

In Fig. [TT1 we plot the variation of the effective grav- 
itational constant G' ff /G e ff versus the redshift z for the 
oscillating scalar field model [ll| with the best fit pa- 
rameters in Table HI In our model, we get G' cS /G c h — 
1.5 x 10~ 13 yr _1 at present. Because of the small value of 
fi^Oj the variation of the effective gravitational constant 
is modest. 
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TABLE I. Results of observational tests from the combination of SN la, CMB, and BAO data. 

Model Best fit parameters x 2 AAIC ABIC 

ACDM ^fi m , = 0.28 466/73 O00 O00 

Oscillating Scalar Field [11] Hl m ,o = 0.28, n , o = 0.01, £ = -25 465.86 3.13 11.11 

a We fix the other parameters at Ho = 72 km s — 1 Mpc _1 . 

b We fix the other parameters at Hq = 72 km s _1 Mpc _1 , m$ = 3.2 X 10~ 31 /i eV, and q = 0. 



6xl0" 13 




0.005 0.01 0.015 0.02 

Z 

FIG. 11. Variation of the effective gravitational constant 
Gc'fl/Geff versus the redshift z for the oscillating scalar field 
model with the best fit parameters in Table[l](The parameters 
same as Fig. I10[) . 

Observationally, there are severe constraints on this 
quantity. From lunar laser ranging tests (3l| . the con- 
straint is reported to be 

^£ = (4 ±9) x 1CT 13 yr- 1 . (33) 

&eff 

Although it is known that the temporal change of the 
gravitational constant affects solar evolution, our predic- 
tion is not inconsistent with the bound from the obser- 
vation. 



IV. CONCLUSION 

According to a Fourier analysis by Hartnett & Hirano 
0, the N—z relation indicates that galaxies have pre- 
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